Probing a ferromagnetic critical regime using nonlinear susceptibility 
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The second order para— ferromagnetic phase transition in a series of amorphous alloys 
(Fe5CosoNii7- x Cr x BieSii2) is investigated using nonlinear susceptibility. A simple molecular 
field treatment for the critical region shows that the third order suceptibility (X3) diverges on both 
sides of the transition temperature, and changes sign at Tc- This critical behaviour is observed 
experimentally in this series of amorphous ferromagnets, and the related assymptotic critical expo- 
nents are calculated. It is shown that using the proper scaling equations, all the exponents necessary 
for a complete characterization of the phase transition can be determined using linear and nonlin- 
ear susceptiblity measurements alone. Using meticulous nonlinear susceptibility measurements, it 
is shown that at times X3 can be more sensitive than the linear susceptibility (xi) in unravelling 
the magnetism of ferromagnetic spin systems. A new technique for accurately determining Tc is 
discussed, which makes use of the functional form of \3 in the critical region. 

PACS numbers: 64.60.Fr, 75.30.Cr, 75.40.Gb, 75.40.Cx 
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I. INTRODUCTION: 

Nonlinear effects become more pronounced in the 
vicinity of a phase transition, and the magnetisation (M) 
measured in the presence of an exciting field (H) can be 
written as 



M = M + X iH + X2H 2 + X3 H 3 



(1) 



Where xi ; \2 and X3 are the first, second and third order 
susceptibilities respectively. The higher order suscepti- 
bilities contain a wealth of information, but owing to the 
fact that they are usually a couple of orders smaller in 
magnitude than the first order susceptibility, are corre- 
spondingly more difficult to measure. 
The second order susceptibility (X2) arises due to the 
presence of a symmetry breaking field and is observed 
only in the presence of a DC field, or in materials with a 
large permanent magnetisation. Though relatively unex- 
plored, X2 has been used to show the coexistence of spin 
glass phase and ferromagnetic order in a reentrant spin 
glass system 

The third order susceptibility (X3) has been commonly 
employed to characterise the spin glass transition, where 
a negative cusp in \3 was used to identify the onset of a 
spin glass transition and was used as a direct probe of the 
divergence of the Edward- Anderson order parameter 0. 
X3 has also been used to distinguish a spin glass from 
a superparamagnetic 0, @- Among long range or- 
der systems, though X3 deduced from dc magnetisation 
measurements has been used for studying the effect of 
quadrupolar interactions in rare earth compounds |(J, it is 
interesting to note that nonlinear susceptibility measure- 
ments on ferromagnetic 0, and antifcrromagnetic0 
materials are not common, and very few reports exist in 
literature. Moreover , most of the previous experimen- 
tal works on nonlinear susceptibility in long range order 
systems are predominantly qualitative and not quantita- 
tive in nature. As far as the study of phase transitions is 



concerned, though the second order paramagnetic to fer- 
romagnetic phase transition has been extensively studied 
using linear susceptibility [lfj; the behavior of X3 across 
the transition has been relatively less explored. 
To bring out the advantages as well as the limitations 
of nonlinear susceptibility measurements in the inves- 
tigation of a para-ferromagnetic phase transition, we 
present a comprehensive study on a set of Fe-Co-Ni-Cr- 
B-Si amorphous alloys. These set of alloys are chosen 
because they were well studied in the past, and it was 
found through linear ac susceptibility and high field dc 
magnetisation measurements that they are reason ably 
well behaved 3-D Heisenberg ferromagnets^J, 0, Il3| . 
As such, amorphous ferromagnets have been extensively 
studied using linear susceptibility to understand the in- 
fluence of quenched disorder on the critical behaviour 
of spin systems. The effect of disorder as seen in the 
values of the critical exponents and the width of the as- 
symptotic critical regime have been studied and reviewed 
extensively 1 14j. Nonlinear susceptibility as measured in 
an amorphous ferromagnet has been reported before [l5|. 
but to the best of our knowledge, this is the first study on 
atleast this class of materials, where linear and nonlinear 
susceptibility alone is used for a complete characterisa- 
tion of the critical regime. 

At the outset, we present a simple molecular field treat- 
ment to show that X3 diverges on both sides of the tran- 
sition, and changes sign at the transition temperature. 
Thus the criticality is observed experimentally, and the 
related asymptotic critical exponents are calculated. It 
has been shown that using the relevant scaling equations, 
all the necessary exponents, required for a complete char- 
acterization of the phase transition can be determined 
using linear and nonlinear susceptibility measurements 
alone. This can be considered as a major advantage of 
measurement of nonlinear susceptibility. Also, owing to 
the functional form of X3 across the phase transition, 
the determination of the true (zero field) transition tem- 
perature (Tc) can be made to a better accuracy. It is 
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shown that subtle features not discernable using normal 
linear susceptibility measurements, can be picked up us- 
ing nonlinear susceptibility, making it a very useful tool 
for studies of spin systems in the critical region. Hence, 
in this paper, we show the following advantages of X3 : 

(i) Determination of all the relevant critical exponents 

(ii) Accurate determination of T c 

(iii) Detecting subtle features not seen in linear suscepti- 
bility measurements. 

It is important to note that though the present study 
deals with amorphous ferromagnets, criticality in X3 can 
be used equally effectively in the study of phase transi- 
tions in their crystalline counterparts. 



ribbon was stuck on to a sapphire single crystal, which 
made up the sample holder, with the temperature sensor 
and a non magnetic heater on the same crystal. This was 
done to avoid temperature gradients across the samples. 
A calibrated Platinum resistance thermometer (PRT) 
was used for measuring the temperature and control to an 
accuracy of 0.01K was achieved using a Lakeshore tem- 
perature controller DRC-93A. Exciting fields from about 
100 mOe to a few Oe and frequencies ranging from 73 Hz 
to 1.33kHz were used for different experimental runs. All 
the measurements were performed with the applied field 
parallel to the long edges of the samples to minimise the 
effect of the demagnetisation fields. 



II. EXPERIMENTAL DETAILS 

A well characterised and extensively studied series of 
amorphous ferromagnets (Fe^Co^NiY!- x Cr x B\^Si\2) 
with x = 5, 10, 15 referred to as A2, A3 and A4 have 
been used in this study. These samples prepared by melt 
quenching are in the form of ribbons of approximate di- 
mensions (10mm X 1mm X 0.03mm). ^lj 
In one of the earlier wor ks 6J, X3 was determined by fit- 
ting the measured isothermal magnetisation using the 
equation 1, or alternatively by plotting M/H vs H 2 , and 
calculating the slope in the linear low field range. 
Another technique involves using the Arrots plots, M 2 
vs H/M , where the slope in the low field range can be 
related to X3i through the relation 

Xs = -Xi\dM 2 /d{H/M)) (2) 

However, these methods are unsuitable for studying the 
critical behaviour across the ferromagnetic transition pri- 
mary due to two reasons; 

(i) The measuring fields are relatively higher, and this 
can at times mask the true critical behaviour of the sys- 
tem 

(ii) X3 is determined after fitting the dc magnetisation 
data. This method is not only time consuming, but also 
introduces significant errors in the value of X3. 

In the present study the linear and nonlinear suscepti- 
bilities are measured by monitoring the change in the 
induced voltage across two oppositely wound secondary 
coils in a home made susceptometer [l6j . Both the linear 
and nonlinear susceptibility was determined by using a 
Lock-in-Amplifier to measure the signal seen at to and 
3 uj respectively, where to is the exciting frequency. Ne- 
glecting the higher order terms, the voltage measured at 
3w is given by 

3 3 

V 3 u> = ^X3# 

where H is the magnitude of the applied field. 

It is normal practice to bundle a few ribbons together to 

increase the measured signal, but in this case, a solitary 



III. RESULTS AND DISCUSSIONS 

A. Linear Susceptibility 

AC susceptibility is amongst the finest tools for the 
study of critical phenomenon, not only due to the rela- 
tive ease with which measurements can be done but also 
due to the fact that very low field measurements can 
be performed. The most important advantage is that it 
directly gives the true initial susceptibility (xo)j which 
otherwise has to be calculated from data taken at com- 
paratively higher measuring fields. Considering that fact 
that the ideal second order paramagnetic-ferromagnetic 
phase transition is defined in zero external field, ac sus- 
ceptibility becomes a very good tool for studying crit- 
icality, as larger exciting fields used during the course 
of measurements are likely to mask the true critical be- 
haviour of the system under study. Also, the Kouvel 
Fisher [K-F] analysis (l7| of the ac susceptibility data is 
a good means of accurately determining both the transi- 
tion temperature (Tq) as well as the susceptibility expo- 
nent (7). Figures la, lb and lc show the Kouvel Fisher 
plots (l/xod/dt(xo^ 1 ) Vs T) for the three samples. The 
inverse of the slopes directly gives us the value of the sus- 
ceptibility exponent 7 and the intercepts of the straight 
line provides us with an independent estimate of the tran- 
sition temperature Tc- This series has been well studied 
in the past|llj|. and the values of Tc and 7 reported are in 
reasonable agreement with the values we have obtained, 
as is shown in Table 1. 



B. Nonlinear Susceptibility 

1 . Theory 

The critical behaviour of X3 nas been investi gate d in 
the past using the Shcrrington-Kirkpatric model|la] and 
the Bethe approximation[19j. The behaviour of X3 wa s 
shown by a simple molecular field theory V by T. Sato 
and Y. Miyako. Here we present a detailed calculation 
based on the molecular field approach. 
Consider a collection of n particles, each with spin 1/2. 
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TABLE I: Values of the susceptibility exponent (7) and the 
transition temperature Tc as determined from the Kouvel- 
Fischer analysis of the first order susceptibility. The values 
determined by the earlier workers is given for the sake of com- 
parison. 





A2 


A3 


A4 




267.44 


222.76 


174.37 




267" 


222. 2 a 


174 6 


7 


1.16 ±0.008 


1.388 ±0.01 


1.41 ±0.01 




1.19 a 


1.38" 


1.73 6 



"Values determined in ref 11] using ac-susceptibility measure- 
ments 

''Values determined in ref ll]using high field dc magnetisation 



When an external field H is applied, the net magnetic 
moment is given by 



m = n,tanh(^) 
kT 



(3) 



where /i is the Bohr's Magneton and k is the Boltzmans 
constant. The net applied field which the dipole sees is 
given by 

H = H + AM 

where A is the molecular field constant. Thus, the mag- 
netisation can now be written as 



M = Nu tanh -^-(H + AM)). 
kT 

where N is the particle density (N = n/V). 



(4) 



using the expression for Tc from the mean field theory, 
kT c = N^ 2 \ 

and comparing the co efficients of H, H 3 , etc we get 



Xi 



1 



X3 = 



T 



3fc 3 T 3 T c ( 



To 



iy 



(7) 



(8) 



It is clear that in the limit T Tc + , Xi diverges 
positively, whereas \3 shows a negative divergence. 

Case B: T < T c 

In the ferromagnetic region, the emergence of spon- 
taneous magnetisation causes M to loose its inversion 
symmetry with respect to H, and now the magnetisation 
can be written as 



M = M + xiH + X2H 2 + X3H 3 



(9) 



where M is the spontaneous magnetisation. Substitut- 
ing this in equation 4 we get 



M = ^ t anh[^;(AM ±(Axi±l)i/±Ax 2 i/ 2 ±A X 3-ff 3 )] 

Expanding tanh(a;) and comparing the coefficients of the 
different powers of H, we get 



Case A: T > T c 

In the paramagnetic region, as the magnetisation (M) 
has inversion symmetry with respect to the applied field 
H, the magnetisation can be written as 



(5) 



M = X iH + X3 H 3 + xsH 5 + .... 
Substituting this in equation 4, we get 

X1H+X3H 3 = Nntanh-^-(H+\lxiH+x 3 H 3 +...}) (6) 
kl 

Now, using the expansion 



tanh(x) = x — —x 3 ; 

o 



M = N^l 



Xi 



X2 



Nu 2 



H\M l/i 3 A 3 M 3 
~~kf 3 k 3 T 3 J 



[W^) 2 (jv^) 2 ] 



kT c 



[£-(l-(^) 2 (^) 2 ) 



kT c M oX i 



N 3 fj, 4 



[1-(^) 2 (^) 2 P 



(11) 



(12) 



(13) 



[1 - (^) 2 (f^) 2 ] 4 N2 V? T (^) 2 [l-(^) 2 (fi) 2 ] 3T C J [ ' 



I 

Equation 11 is solved for the spontaneous magnetisation 

and yields for T < Tc Substituting this value of (If 2 -) 2 in the value of Xl an d 



( ^)2 = 3( ^ ) 2 (1 _^L ) (15) 

N/i y T c T c 
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simplifying, we get 



Xi 



jy [3^ - 2] 
MTc [1 - 



(16) 



Hence it is clear that xi nas a positive divergence as 

Now, substituting equations 15 and 16 in the value of %3, 
and simplifying, we obtain 



X3 



(JLtL\± 



AT 



3(1-^)4 {NnfTc 



(17) 



where it is obvious that like \i ? X3 too diverges positively 
in the limit T => T c~ 



2. Experiments 

Normally, ^3 is seen to show a sharp negative peak 
around the transition. However, as is clear from the pre- 
vious section, theroretically, X3 is expected to diverge on 
both sides of the phase transion, and change sign exactly 
at Tq- However, this critical behaviour is not easy to 
observe. The reason for it is that the expansion of M in 
terms of higher powers of H, which was shown in the pre- 
vious section is truely speaking only valid for small values 
of H. Hence, to observe the true critical behaviour, mea- 
surements have to be done at very low measuring fields, 
and higher fields tend to smear off the transition. It is 
difficult to predict apriori the fields at which this criti- 
cal behaviour will be seen, as it depends on the extent 
of nonlinearity in the system, and thus differs from one 
sample to another. It must be kept in mind that at higher 
measuring fields, the contribution from domains can be- 
come large, thus masking the true critical behaviour of 
the spin system. 

Careful low field measurements have shown this critical 
feature in this series of samples. Fig 2a, 2b and 2c shows 
the critical behaviour in the samples A2, A3 and A4 
respectively. As predicted in the previous section, \3 
changes sign across Tc and the exact Tc can be directly 
determined from the crossover of X3 m the temperature 
axis. This experimental Tc values are given in Table 2. 
It is to be noted here that the Tc determined from the 
crossover of X3 matches well with that determined from 
the Kouvel-Fischer plot (Table 1). This matching of Tc's 
from two different measurements not only substantiate 
the fact that the measured X3 is a genuine response of 
the spin system, but also that the crossover of X3 can be 
used as a direct method to determine the Tc- The other 
advantages of determination of Tc from the crossover of 
X3 will be discussed in the next section. 
As is seen in fig 2a, a double transition showing two 
crossovers is observed in the sample A2. As will be shown 
later, this sample gives an uncharacteristically low value 
of the susceptibility exponent^), a fact observed by the 



earlier workers as well. Considering the fact that this 
sample has the largest percentage of Ni and is closest to 
the critical concentration^), this double transition is 
not difficult to explain, as in this class of materials, the 
formation of clusters showing a distribution of Tc's have 
been speculated^J0]. However, it is interesting to note 
that no direct evidence of such a distribution is evident 
from the linear susceptibility measurements. This only 
goes on to show that low field nonlinear susceptibility is a 
more sensitive tool than the linear susceptibility, as far as 
studies of spin systems near the transition is concerned. 
As mentioned earlier, the field in which the measurement 
is performed, is an important consideration in the study 
of criticality in X3- This is clearly shown in fig 3, which 
shows the field dependence of X3 i n the sample A3 at a 
measuring frequency of 133.33 Hz. It is clearly seen that 
the critical feature is sharpest at the lowest measured 
field, and increasing the measuring field wipes off this 
critical behaviour. It is interesting to note that we have 
observed a non trivial frequency dependence in this series 
of samples . This could possibly be due to the presence of 
superparamagnetic clusters, as no frequency dependence 
would be expected in a long range ferromagnetically or- 
dered system. These results will be dealt with in a later 
communication. 

Large uncertainties in ^(X3 _1 ) make it very difficult to 
determine the value of the critical exponent associated 
with the third order susceptibility (73), using the Kouvel- 
Fischer formalism. Hence the value of 73 is determined 
by plotting a double logarithmic plot of (— |)x3-ff 2 Vs e 
where e = ( T ~ Tc ) is calculated by taking the crossover 
point of X3 as Tc- Fig 4a, 4b and 4c shows the double 
logarithmic plots for samples A2, A3 and A4 respectively. 
As shown in the figure, as T — » Tc + , a curvature from 
the straight line fit is clearly seen. This arises due to 
the presence of higher order terms close to the transition 
temperature. Exponent calculations gives the value of 73 
to be 4.57, 4.80 and 5.04, for this series of samples, which 
matches well with the 3-D Heisenberg value of 4.88. 
Using the scaling relations 
73+ = 7 + 2A ; 

A = 7 + A 

S = 1 + 7//8; and 

a + 2/3 + 7 = 2, 

(where a ,/3 ,7 and A are the exponents associated 
with the specific heat, magnetisation, susceptibility and 
the gap exponent [23] respectively) all the exponents, re- 
quired for a complete characterisation of the system have 
been determined, and is given in table 2. As is clearly 
seen our values match reasonably with that given in 
ref. 11] The only exponent in which a large difference 
is seen is a which is calculated by using the Rushbrooke 
equality, (a + 2(3 + 7 = 2). The ambiguity in the value of 
this exponent is probably because the values determined 
in ref[llj is calculated by using both high dc field and low 
ac field measurements, as is the general practice, whereas 
in our case a is determined by using low field linear and 
nonlinear ac susceptibility measurements alone. 
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TABLE II: Values of the all the critical exponents as calcu- 
lated using linear and nonlinear ac susceptibility in this series 
of samples. Tc(xs) is the value of the transition temperature 
determined from the crossover of Xd, m the temperature axis. 
The values quoted by the previous workers is also given for 
the sake of comparison 





A2 


A3 


A4 


7 


1.16 ±0.008 


1.388 ±0.01 


1.41 ±0.01 




1.19" 


1.38" 




73 


4.57 ±0.08 


4.805 ± 0.12 


5.04 ± 0.09 


A 


1.705 


1.708 


1.815 


P 


0.545 


0.32 


0.405 




0.35 a 


0.41 a 


0.52 a 


a 


-0.25 


0.03 


-0.22 




0.2 a 


-0.2 a 


-0.7 a 


S 


3.12 


5.33 


4.48 




4.42 a 


4.49 a 


4.32 a 


Tc(x3) 


267.05 


222.55 


173.45 



"From Ref £3 



C. Determination of Tc using Nonlinear 
susceptibility 

Accurate determination of Tq has always been an im- 
portant consideration in the study of criticality across the 
ferromagnetic transition. The most popular techniques 
used till date have been the K-F formalism (discussed in 
the earlier section) and the Arrots plots plj. 
Arrots plots works on the principle of plotting the isother- 
mal experimental data as M 2 Vs H/M at different 
temperatures in the transition region. They should be 
straight lines in the critical region with the intercepts of 
these lines with the (H/M) axis being positive if T > Tc, 
and negative when T < Tc- However, it is important 
to note that this is valid only when the domain allign- 
ment is complete, ie the low field experimental data is 
completely excluded. This can at times be a non trivial 
problem, specially in systems with large field induced ef- 
fects, causing uncertainties in the accurate determination 
of Tc- Also, in inhomogeneous systems, (like amorphous 
ferromagnets), the Arrots plots, tend to be curved, even 
at large measuring fields, thus making the identification 
of Tc even more difficult . 

As has been described earlier, \3 diverges on both sides 
of the transition temperature, and changes sign exactly 
at Tc- Moreover, this critical behaviour is extremely 
sensitive to the magnitude of the field in which the mea- 
surement is done, and this critical feature is seen to van- 
ish with increasing field. This makes X3 a very promising 



candidate for an accurate determination of Tc, as even 
by theory, the ferromagnetic transition is defined for the 
limit H =>• 0, and any large measurement field is likely 
to change the thermodynamic properties of the system 
under study. Plotting the temperature of crossover (Tc) 
Vs the field used during the measurement (H ac ), it is 
clearly seen that the Tc varies linearly with the applied 
field. Thus, extrapolating this to H = 0, one will be 
able to determine Tc to a very high accuracy. This is 
shown in figure 5, where Tc vs H is plotted for the sam- 
ple A3. The intercept at H = 0, gives us the value of 
T c = 222.43K. 



IV. CONCLUSIONS 

A simple molecular treatment of the second order para- 
ferromagnetic phase transition shows that X3 diverges on 
both sides of the transition, and changes over sign at 
Tc- This criticality has been experimentally observed in 
a well studied series of amorphous alloys, and the related 
assymptotic critical exponents have been determined. It 
is shown that sensitive nonlinear susceptibility measure- 
ments , along with linear susceptibility, can be used to 
fully characterize the phase transition, without the need 
for any other type of measurements, like dc magnetisa- 
tion, specific heat, etc. Moreover, due to the functional 
form of %3 at the transition, the transition temperature 
Tc can be determined accurately. Using the field de- 
pendence of X3, an extrapolation to true zero field can 
be made, which gives us a unique way of determining 
the transition temperature. Subtle features not seen in 
the first order susceptibility, like the presence of magnetic 
clusters with a distribution of Tc's is observed. This only 
goes on to show that at times nonlinear susceptibility can 
be a more sensitive tool than their linear counterparts in 
understanding the magnetism of spin systems in the crit- 
ical region. 
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FIG. 1: Kouvel Fischer Plots of the real part of the first or- 
der susceptibility for the samples A2, A3 and A4 respectively. 
The inverse of the slope gives the value of the susceptibil- 
ity exponent (7), and the intercept on the temperature axis 
provides the transition temperature (Tc)- 

FIG. 2: The critical feature of Xi as observed in the real part 
of the third order susceptibility in the samples A2, A3 and 
A4 respectively. The measurements were done at a frequency 
of 133.33 Hz and the exciting field is of the order of 200 mOc. 

FIG. 3: The field dependence of the third order susceptibility 
as seen in the sample A3 at a measuring frequency of 133.33 
Hz. As is clearly seen, the critical behaviour is supressed with 
increasing applied field. 

FIG. 4: Double logarithmic plots of \3 versus the reduced 
temperature for the samples A2, A3 and A4 respectively. The 
slopes of the straight lines give the value of the exponent 73. 

FIG. 5: The temperature of crossover Tc versus the mea- 
suring field Hac for the sample A3. A linear behaviour is 
observed, and the interpolation to zero applied field gives us 
a unique way of determining the true transition temperature. 
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